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Abstract 

In conventional approaches to the homogenization of random particulate composites, the component 
phase particles are often treated mathematically as vanishingly small, point -like entities. The elec- 
tromagnetic responses of these component phase particles are provided by depolarization dyadics 
which derive from the singularity of the corresponding dyadic Green functions. Through neglecting 
the spatial extent of the depolarization region, important information may be lost, particularly 
relating to coherent scattering losses. We present an extension to the strong-property-fluctuation 
theory in which depolarization regions of nonzero volume and ellipsoidal geometry are accommo- 
dated. Therein, both the size and spatial distribution of the component phase particles are taken 
into account. The analysis is developed within the most general linear setting of bianisotropic 
homogenized composite mediums (HCMs). Numerical studies of the constitutive parameters are 
presented for representative examples of HCM; both Lorentz-reciprocal and Lorentz-nonreciprocal 
HCMs are considered. These studies reveal that estimates of the HCM constitutive parameters in 
relation to volume fraction, particle eccentricity, particle orientation and correlation length are all 
significantly influenced by the size of the component phase particles. 
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1 Introduction 

The homogenization of a composite medium comprising two (or more) component phases provides 
the backdrop for this study. The composite may be regarded as an effectively homogenous medium 
as long as wavelengths are sufficiently long compared with the dimensions of the component phase 
particles. In electromagnetics, the estimation of the constitutive parameters of homogenized com- 
posite mediums (HCMs) is a matter of long-standing, and ongoing, scientific and technological 
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importance [1]. Indeed, recent developments pertaining to HCM-based metamaterials serve to 
highlight the need for accurate formalisms in order to estimate the constitutive parameters of 
complex HCMs [2, 3]. 

In conventional approaches to homogenization the microstructural details of the component phases 
are often inadequately incorporated [4]. For example, the standard versions of the widely-applied 
Maxwell Garnett [5] and Bruggeman [6] formalisms utilize simplistic descriptions of the spatial dis- 
tributions and sizes of the component phase particles. An alternative approach to homogenization is 
provided by the strong-property-fluctuation theory (SPFT), in which a comprehensive description 
of the distributional statistics of the component phases can be accommodated. While the origins 
of the SPFT lie in wave propagation studies for continuous random mediums [7, 8], the SPFT has 
lately gained prominence in the homogenization of particulate composites [9, 10, 11, 12, 13]. By 
means of the SPFT, the HCM constitutive parameters are estimated as successive refinements to 
the constitutive parameters of a homogenous comparison medium. Iterates are expressed in terms 
of correlation functions describing the spatial distributions of the component phases. In principle, 
correlation functions of arbitrarily high order may be incorporated; in practice, the SPFT is usually 
implemented at the second order level of approximation. In fact, convergence of the SPFT scheme 
at the second order level of approximation has been established for a wide range of linear HCMs 
[14]. Within the second order SPFT, a two-point correlation function and its associated correla- 
tion length characterize the component phase distributions. At lowest order (i.e., zeroth and first 
order), the SPFT estimate of HCM constitutive parameters is the same as that of the Bruggeman 
homogenization formalism [9, 13]. 

The size of the component phase particles may be explicitly incorporated within homogenization 
formalisms via depolarization dyadics [4]. These dyadics are central to homogenization analyses 
as they characterize the electromagnetic field inside component phase particles embedded within 
a homogenous background. Often the component phase particles are treated as vanishingly small, 
point-like entities. Thereby, the corresponding depolarization dyadic is represented by the singular- 
ity of the associated dyadic Green function [15, 16]. However, potentially important information is 
lost through neglecting the spatial extent of the component phase particles — especially if coherent 
scattering losses are under consideration [17, 18]. To address this issue, extended versions of both 
the Maxwell Garnett formalism [19, 20, 21] and the Bruggeman formalism [21, 22] have been de- 
veloped in which a nonzero volume is attributed to the component phase particles. However, these 
analyses apply only to isotropic HCMs and adopt a simplistic description of the distributional 
statistics of the component phases. Recently, an extended version of the SPFT was established 
for anisotropic dielectric HCMs, in which the sizes of the component particles and their spatial 
distributions were taken into account [23]. 

In the following sections, the second order SPFT for the most general class of linear HCMs, namely 
bianistropic HCMs [24], is extended to accommodate component phase particles of nonzero size. 
The depolarization dyadic appropriate to an ellipsoidal particle of nonzero size, embedded in a 
bianisotropic medium, is developed in §2. The incorporation of this depolarization dyadic within 
the SPFT is then outlined in §3. The influence of the size of the component phase particles 
upon the estimates of the HCM constitutive parameters is explored in numerical studies in §4, for 
both Lorentz-reciprocal and Lorentz-nonreciprocal HCMs. Lastly, a few concluding remarks are 
provided in §5. 
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The following notation is adopted: Vector quantities are underlined. Double underlining and 
normal (bold) face signifies a 3x3 (6x6) dyadic. The inverse, adjoint, transpose and determinant 
of a dyadic M are denoted by MT 1 , adj [M] , M T and det [M] • The 3x3 (6x6) identity dyadic is 
represented by 7 ( I ) . All field-related quantities are implicitly functions of the angular frequency 
u. The permittivity and permeability of free space are denoted as e and fj, , respectively; the 
free-space wavenumber is k = u;^/e 0/ u„ ; and A = 27r/k . The real and imaginary parts of z G C are 
represented by Re z and Im z, respectively. A compact representation of the constitutive parameters 
for the homogeneous bianisotropic medium specified by the Tellegen constitutive relations [24] 



D{r) = g • E(r) + £ • H(r) 
B(r) = ( ' E{r) +^ • H(r) 



(1) 



is provided by the 6x6 constitutive dyadic 



K 



(2) 



Herein, e and /U are the 3x3 permittivity and permeability dyadics , respectively, while £ and £ 
are the 3x3 magnetoelectric dyadics. In the following, subscripts on K are used to identify the 
particular medium that K describes. 



2 Depolarization region 

Let us consider an ellipsoidal particle of volume Vj 1 , oriented arbitrarily in R 3 . The ellipsoidal 
surface of is parameterized by 

r e (9,<f>)= V g>r(9,4>), (3) 

where r (9, (f>) is the radial unit vector specified by the spherical polar coordinates 9 and (j). The 
3x3 shape dyadic U_, which is real symmetric with unit determinant, maps the spherical region 
V 7 * of radius rj onto the ellipsoidal region V<P. The linear dimensions of the ellipsoidal particle, as 
determined by i], are assumed to be sufficiently small that the electromagnetic long-wavelength 
regime pertains, but not vanishingly small. 

Suppose now that the ellipsoidal particle is embedded within a bianisotropic comparison medium, 
characterized by the 6x6 constitutive dyadic K . The comparison medium is homogeneous. The 
electromagnetic response of the ellipsoidal particle is provided by the depolarization dyadic [16] 

D(r?)= f G (r)d 3 r=[ G (U • r) d 3 r . (4) 

Herein, G_ cm (r) is the 6x6 dyadic Green function of the comparison medium which satisfies the 
nonhomogenous vector Helmholtz equation [16] 

•G (£-/)= I6(r-r!), (5) 



L(V) +iuK 
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with the linear differential operator 



L(V) = 

and S(r — r') being the Dirac delta function. 



g V x £ 
-V x / 



(6) 



Explicit representations of Green functions are not generally available for anisotropic and bian 
isotropic mediums [25]. However, it suffices for our present purposes to consider the Fourier trans- 
form of G (r), namely 

J3 



Jr_ 

which is delivered from equation (5) as 



G (q) = — A (q) 

= cm v — ' [_= cm v — ' 



where 



A (q) = 

= cm v — ' 

Thereby, equation (4) yields [15, 16] 



g {q/uj) x X 

■(?/w)xI 



+ K 



m ^/^(^-co S(W )) icJrl .,)A. 



(7) 

(8) 
(9) 

(10) 



In order to consider the depolarization dydic of an particle of nonzero volume, we express D.(rj) as 
the sum 

D(r ? )=D>°(r ? ) + D°, (11) 
The two terms on the left side of (11) are given by 



with 



G°° (J/" 1 • §) = lim G OJ~ l • q) 

^ij,(0,4>)qq -^{0,4>)qq 



iuj b(0, i 



-a^,4>)qq a e {9,(j))qq 



wherein the scalars 



a p (6,4>) =q -IT 1 -p 'Ur l 'q, 



(P = e,C,£,Ai) 



(12) 
(13) 



(14) 
(15) 



(16) 
(17) 
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and 

b{9, 4>) = [a e (0, <j>) ^(6, <j>)] - [as(6, <j>) a c (0, </>)] . (18) 
The volume integral (13) simplifies to the 77-independent surface integral [15, 16] 

D° = — U- 1 - ( [ I G°° OJ- 1 - q) sinO d9 deb] • U~K (19) 
= 4vr= V^=o Je=o= cm = ~ J = 

For certain Lorentz-reciprocal comparison mediums, the volume integral (12) which yields D >0 (r?) 
may be reduced to a surface integral, but for a general bianisotropic comparison medium no such 
simplifications are available. The integrals (12) and (19) may be evaluated using standard numerical 
techniques [26]. 

The dyadic D° represents the depolarization contribution arising from the vanishingly small region 
of volume limV^, whereas the dyadic D >0 (ry) provides the depolarization contribution arising from 

?7 — >0 — 

the region of nonzero volume I Ve — limV^ 7 J. In homogenization studies, it is common practice 

to neg lect D >0 (r/) and assume that the depolarization dyadic is given by D° alone [4]. However, 
studies of isotropic [17, 18, 19, 21, 22] and anisotropic [23] HCMs have emphasized the importance 
of the nonzero spatial extent of depolarization regions. 



3 Homogenization 

The SPFT may be implemented to estimate the constitutive parameters of HCMs [9]. Let us 
consider the homogenization of a two-phase composite wherein the two component phases, labelled 
as a and b, comprise ellipsoidal particles of shape specified by U and linear dimensions specified 
by rj > 0. A random distribution of identically oriented particles is envisaged, such that all space 
V is partitioned into parts V a and V& containing the phases labelled a and b, respectively. We 
consider the most general linear scenario wherein the component phases a and b are taken to be 
bianisotropic mediums with 6x6 constitutive dyadics K q and K ft , respectively. 

The distributional statistics of the component phases are described in terms of moments of the 
characteristic functions [13] 

f 1, L€V t , 

$i(D = { (£ = a,b). (20) 

I 0, r V e , 

The volume fraction of phase £, namely fa , is given by the first statistical moment of ; i.e., 
(&e(z) ) = fe ■ Clearly, f a + fb = 1- The second statistical moment of $v provides a two-point 
covariance function. We adopt the physically-motivated form [27] 



(Ml) Ml')) 



( (MD)(Ml')), \UT l • {L-r')\ >L, 

(21) 

I (Ml)), \g~ l • (l-l')\<l, 



where L > is the correlation length, which is taken to be much smaller than the electromagnetic 
wavelengths. Over a range of physically-plausible covariance functions, it has been shown that the 
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specific form of the covariance function has only a secondary influence on SPFT estimates of HCM 
constitutive parameters [28]. 



3.1 Zeroth order SPFT 

The nth order SPFT estimate of the HCM constitutive dyadic, namely K# CM > is based upon the 
iterative refinement of the comparison medium constitutive dyadic, namely K cm - To zeroth order 
and first order, the SPFT permittivity estimate is identical to the comparison medium permittivity 
[13]; i.e., 

K'l M = KW „ = K . (22) 

= HCM =HCM =cm y ' 

The well-known Bruggeman homogenization formalism provides the estimate of K cm [13]. That 
is, K cm emerges through solving the nonlinear equations 

faX a/Cm + fbK b/Cm =^ ( 23 ) 

wherein the polarizability density dyadics 

t jCm = --(I*-*™) • [1 + ^U(V)' (K,-Kj]"\ (£ = a,b). (24) 



3.2 Second order SPFT 



The SPFT is most widely implemented at the second order level (also known as the bilocal approx- 
imation) which provides the following estimate of the HCM constitutive dyadic [13] 



= HCM =cm 



K 



iuj 



I + SPl .D(r/)1 -£ [2] . 



(25) 



Thus, the particle size rj influences K'| CM directly through the depolarization dyadic D(r?) and 
indirectly through the mass operator [7] dyadic term 

I [2] = fah (V /cm - X b/cm ) ' D >0 (L) • (r' cm - ^ cm ) . (26) 

Notice that the correlation length L — which plays a key role in the second order SPFT — does 
not feature in the zeroth order SPFT. 



4 Numerical studies 



We now apply the theoretical results presented in Sections 2 and 3 to two specific bianisotropic 
homogenizations scenarios: in §4.1 a biaxial bianisotropic HCM is considered and in §4.2 a Faraday 
chiral medium [29, 30] is considered. The HCM in §4.1 is Lorentz-reciprocal [31] whereas the 
HCM in §4.2 is not. Numerical studies are presented for representative examples, in order to 
explore the influence of rj in relation to volume fraction, particle eccentricity, particle orientation 
and correlation length. In view of the vast parameter space associated with bianisotropic mediums, 
only an illustrative selection of graphical results are provided here. 
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The following calculations were carried using an angular frequency u = 2tt x 10 10 rads 1 . Hence, 
A = 2ir/k = 0.030 m. 



4.1 Biaxial bianisotropic HCM 

The homogenization of (i) a biaxial dielectric medium described by the constitutive dyadic 

" e diag (e*, e y a , e z a ) g 
= a ~ [ g fi {> l 

and (ii) an isotropic chiral medium described by the constitutive dyadic 

e e b l V e oMo6I 



(27) 



-iy/eoHo&L Mo I 



(28) 



is investigated. The constitutive parameters selected for calculations are: t x a = 4 + i0.12, e« 
3 + iO.l, = 1.5 + i0.08; e 6 = 2.5 + iO.l, 6 = 1 + *0.07 and p b = 1.75 + z0.09. 



The shape dyadic of the constituent particles is taken to be 

1 



U 



Riip) • [diag(C4,^,^)] 



with 



cos ip sin (/? 
R r (tp) = —sin 99 cos if 
1 



(29) 



(30) 



Thus, the principal axes of the ellipsoidal particles lie in the xy plane rotated by an angle 92, and 
along the z axis. The shape parameters selected for calculations are: U x = 1 + p, U y = 1 and 
U z = 1 — 0.5/5; the eccentricity of the ellipsoids is varied through the parameter p. 

The corresponding HCM is a Lorentz-reciprocal, biaxial, bianisotropic medium. In this case, the 
volume integral in (12) for D >0 (??) is analogous to one which arises in the development of the second 
order SPFT [13]. Thus, we may express D >0 (f?) as the surface integral 



,.,3 ft 1 r p lr !1 r 
— 8tti 7 0=o J e=0 - K_ l ^ L 

+N(-£T 1 * 9)1 1""^ + _?_N(0) I sin^ dO #, 

mur 1 -3) = K^j h i^jur 1 • z>] - [k c jg- 1 • <o\ cjr 1 • §)} 

and k± being the roots of det ^^(C/ -1 • q) 



with 



(31) 



(32) 
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4.1.1 Particle size and volume fraction 



To focus upon the effect of particle size 77 in relation to volume fraction f a , we set the eccentricity 
p = 0, the orientation angle ip = and the correlation length L = 0. The corresponding HCM 
constitutive dyadic has the form 

e diag (e HCM , ^ucm-> € hcm) ^v 76 ' 1 ^ {%hcm-> ^hcm^hcm) 
—i^/e p diag (£# CM , ^hcm^hcm) Mo diag (p, HCM , ^hcm^hcm) 



= hcm 



In Figure 1, the real and imaginary parts of € x ^q Z m are plotted against / a G (0,1) for r//A € 
{0, 0.05, 0.1}. Notice that the HCM parameters are constrained to coincide with those of component 
phase b in the limit f a — ► 0, and those of component phase a in the limit f a — > 1. The influence of 
77 is more obviously observed on the imaginary parts of t^Q M than on the real parts. Indeed, the 
imaginary parts of e^, z M for 77/ A = 0.1 at mid-range values of / a are approximately twice as large 
as they are for rj/X = 0. The corresponding graphs for the HCM magnetoelectric parameters £#c»m 
and permeability parameters P^cm are qualitatively similar to those presented for the permeability 
parameters in Figure 1. 



4.1.2 Particle size and particle eccentricity 

Next we turn to the effect of particle size 77 in relation to particle eccentricity, as specified by p. 
The volume fraction is fixed at f a = 0.5, the orientation angle at ip = and the correlation length 
at L = 0. The corresponding HCM constitutive dyadic has the form (33). 

As a typical example of the behaviour of HCM constitutive parameters, the real and imaginary 
parts of ^hcm are graphed versus p G (0, 1) for r]/\ G {0, 0.05, 0.1} in Figure 2. Regardless of the 
value of 77, the constitutive parameters vary substantially — particularly their imaginary parts — as 
p increases. Furthermore, there are significant differences in the plots of ^}jcm presented for the 
three values of 77. The most striking differences are observed in the plots of the imaginary parts of 
£hcm- The corresponding graphs for the constitutive parameters not presented in Figure 2 (i.e., 
e HCM' €hcm an d P^hcm) are broadly similar to those given in Figure 2. 



4.1.3 Particle size and particle orientation 



In order to investigate the effect of particle size 77 in relation to particle orientation, we fix the 
volume fraction f a = 0.5, particle eccentricity p = 1 and the correlation length L = 0. The 
resulting HCM has a constitutive dyadic of the form 



= HCM 



I 



€ HCM 
HCM 





e HCM 

y 

€ HCM 




HCM 
%HCM 





^HCM 



\ 


e HCM J 






^hcm 

£hcm J 



Mo 




• (34) 
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Illustrative numerical results are displayed in Figure 3, wherein the real and imaginary parts of 
I^HCM are Plotted against (p G (0,7r/2) for r)/\ £ {0,0.05,0.1}. The off-diagonal constitutive 
parameter // vanishes in the limits (p — ► and 7r/2 (as do e* and £*). Both the real and imaginary 
parts of /x* are strongly influenced by the particle size 77, especially for mid-range values of </?. The 
diagonal constitutive parameter \i v is also clearly sensitive to rj. In the case of fi y , the differences 
in behaviour for the three values of rj are most apparent as ip approaches and ir/2. The graphs of 
//* are symmetric about cp = 7r/4, but those of fi y are not. The HCM constitutive parameters that 
are not represented in Figure 3 (i.e., e#cM> £hcm anc ^ V^hcm ) exhibit behaviour with respect to 
93 which is generally similar to that exhibited by m Fig ure 3. 



4.1.4 Particle size and correlation length 

Lastly in this section, particle size rj is considered in relation to correlation length L. To do so, 
the following parameters are fixed: volume fraction f a = 0.5, orientation angle ip = and the 
eccentricity p = 0. The constitutive dyadic of the HCM which arises has the form (33). 

In Figure 4, graphs of the real and imaginary parts of e x HCM versus k L G (0, 0.2) are provided for 
rj/L G {0,0.5,0.95}. It is clear that the imaginary part of e'jj CM is strongly affected by increasing 
L; the real part of Cjjcm ^ s a ^ so affected but to a lesser degree. Furthermore, Cjjcm * s mucri more 
sensitive to L at larger values of 77. The behaviour observed in Figure 4 for e x HCM with respect to 
L is also generally observed in the HCM constitutive parameters £jf CM , £hcm ano - I^hcm wn ich 
are not represented in Figure 4. 



4.2 Faraday chiral medium 



For our second homogenization scenario, we explore the homogenization of (i) a gyrotropic magnetic 
medium described by the constitutive dyadic 



K 



e e a L 







V 



Ma 
■ 





■ 9 




\ 


Ma / 



(35) 



and (ii) an isotropic chiral medium described by the constitutive dyadic (28). The constitutive 
parameters selected for calculations are: e a = 1.2 + i0.02, = 3.5 + i0.08, yf a = 1.8 + i0.05, 
= 1.4 + i0.04; e b = 2.5 + iO.l, £ 6 = 1 + i0.07 and // 6 = 1.75 + i0.09. As in §4.1, the shape dyadic 
of the constituent particles is taken to have the form (29), with the shape parameters selected for 
calculations being: U x = 1 + p, U y = 1 and U z = 1 — 0.5p. 

The HCM which results is a Faraday chiral medium [29, 30, 32]. A HCM of the same form also 
arises from the homogenization of a magnetically-biased plasma and an isotropic chiral medium 
[33]. 
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4.2.1 Particle size and volume fraction 



We begin by considering the effect of particle size rj in relation to volume fraction f a . Accordingly, 
the eccentricity is fixed at p = 0, the orientation angle at ip = and the correlation length at L = 0. 
The HCM constitutive dyadic has the form 



= HCM 



( 



^HCM 
~ le HCM 



HCM 

e HCM 




■>HCM 







-HCM 



it? 



HCM 
?~x 

'HCM 






^HCM / 



Mo 



CI 



*HCM 

~^HCM 




3 

HCM 



Phcm 

■ 9 

~ l PHCM 



da; 

^HCM 





l t x HCM 
PhCM 



\ 


£jfCM / 




Phcm 



(36) 



In Figure 5, the real and imaginary parts of Pjj CM are plotted against f a G (0,1) for rj/A G 
{0,0.05,0.1}. As is the case in Figure 1, the HCM constitutive parameters are constrained such 
that they coincide with those of component phase b and a in the limits f a ^0 and 1, respectively. 
The effect of rj on the real parts of p X jf CM are relatively modest. In contrast, rj has a profound 
effect on the imaginary parts of p^CMi especially for mid-range values of f a . The pattern of 
behaviour presented in Figure 5 for the HCM permeability parameters p X ifcM 1S m i rr ored by the 
HCM permittivity parameters ^ c g M and magnetoelectric parameters ^^M' as weu as Phcmi 
which are not displayed in Figure 5. 



4.2.2 Particle size and particle eccentricity 

The effect of particle size r] in relation to particle eccentricity is considered next. We set the volume 
fraction f a = 0.5, orientation angle <p = and the correlation length L = 0. The HCM constitutive 
dyadic then has the form 



= HCM 



€ HCM 



If 



-ie 



9 

HCM 




ex 
HCM 

l ^HCM 





HCM 

y 

C -HCM 



^HCM 
^>HCM 




Po 




which is rather more general than the form (36). 

As a representative example, graphs of the real and imaginary parts of versus P £ (0, 1) are 

exhibited in Figure 6 for r//A G {0,0.05,0.1}. The real part of CfjcM ^ s a strong function of p, 
whereas the imaginary part varies less as p increases. The sensitivity of both the real and imaginary 
parts of Chcm *° P ^ s clearly influenced by the value of r\. Patterns of behaviour with respect to 
p that are qualitatively similar to those displayed in Figure 6 are found for the HCM constitutive 



x,y,z,g f-x,y,gl,g2 



parameters £ HCM ^, £' 



HCM 



and M_ffCM w hidi are not represented in Figure 6. 
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4.2.3 Particle size and particle orientation 



Now we turn to effect of particle size rj in relation to particle orientation. Let us fix the following 
parameters: volume fraction f a = 0.5, particle eccentricity p = 1 and the correlation length L = 0. 
Consequently, the HCM constitutive dyadic has the form 



= HCM 



( 



—le 



HCM 
92 



if 



HCM 





91 

HCM 

y 

e HCM 










'HCM 



£>HCM 

lt ^HCM 




■tf 1 

lt ^HCM 
^HCM 



Chcm 
-iC 92 

^>HCM 





*»1 



^HCM 










Mo 



^HCM 
_ ■ 92 
l ^HCM 



■ 9l 

l l x HCM 

y 

^HCM 



\ 






Vhcm 



(38) 

which is more general than (36) and (37). As illustrative examples, the the real and imaginary parts 
of ^hcm are plotted in Figure 7 against ip € (0,7r/2) for r//A G {0,0.05,0.1}. The particle size rj 
has a strong influence on the real and imaginary parts of e 9 ^ CM , as well as on the imaginary part of 
€ hcm- T ne influence on the real part of e x HCM is smaller by comparison, but still significant. The 
graphs for e 9 ^ CM are symmetric about the line <p = 7r/4 whereas those for e| Ctf are not. Broadly 
similar behaviour is exhibited by the HCM constitutive parameters not plotted in Figure 7, namely, 

,V,z,g2 j-x,y,z,gl,g2 j-x,y,gl,g2 , x,y,z,gl,g2 
HCM' >>HCM ' ^HCM dIlu t^HCM 



4.2.4 Particle size and correlation length 

Finally, to focus upon the effect of particle size rj in relation to correlation length L, the volume 
fraction is fixed at f a = 0.5, the orientation angle at if = and the eccentricity parameter at p = 0. 
The corresponding HCM constitutive dyadic has the form (36). 

In Figure 8, the real and imaginary parts of p g HCM are plotted against k L 6 (0,0.2) for rj/L G 
{0,0.5,0.95}. The pattern of behaviour with respect to L is similar to that presented in Figure 4 
for the biaxial binisotropic HCM. That is, the imaginary part of M/^cm ^ s more obviously sensitive 
to inceasing L than is the real part. In addition, both the real and imaginary parts of p- g HCM are 
more sensitive to L at larger values of 77. The other HCM constitutive parameters, namely ^q 9 m , 
£hcm ana - ^HCM' res P on d in a generally similar manner as L increases for the three values of 77 
considered here. 



5 Concluding remarks 

Homogenization formalisms, such as the widely-used Maxwell Garnett and Bruggeman formalisms, 
often inadequately take into account the distributional statistics and sizes of the component phase 
particles. The SPFT — through describing the distributional statistics of the component phases 
in terms of a hierarchy of spatial correlation functions — provides a conspicuous exception. In 
the preceding sections, an extension to the SPFT for the most general linear class of HCM is 
developed, in which a nonzero volume is attributed to the component phase particles. By means of 
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extensive numerical calculations, based on Lorentz-reciprocal and Lorentz-nonreciprocal HCMs, it 
is demonstrated that estimates of the HCM constitutive parameters in relation to volume fraction, 
particle eccentricity, particle orientation and correlation length are all significantly influenced by 
the size of the component phase particles. 

It is particularly noteworthy that the influence of the particle size is generally stronger on the 
imaginary parts of the HCM constitutive parameters than it is on the corresponding real parts. In 
this respect, the effect of r\ is reminiscent of the effect of the correlation length in the second order 
SPFT [34]. Increasing the correlation length for the second order SPFT generally results in an 
increase in the degree of dissipation associated with the HCM. This dissipative loss is attributed 
to radiative scattering losses from the macroscopic coherent field [13, 35]. It may be observed in 
Figures 4 and 8 (and in numerical results for other HCM constitutive parameters not presented 
here) that the effects of particle size and correlation length on the estimates of the imaginary 
parts of the HCM constitutive parameters are generally cumulative. This suggests that coherent 
scattering losses associated with the HCM become greater as the particle size increases. A similar 
finding was reported for anisotropic dielectric HCMs [23]. 

In conclusion, the importance of incorporating microstructural details, such component particle size 
and spatial distribution, within homogenization formalisms is further emphasized by this study. 
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Figure 1: Real (left) and imaginary (right) parts of the HCM constitutive parameters plotted 
against volume fraction f a 6 (0, 1) for rj/\ = (dashed curves), f]/X = 0.05 (broken dashed curves) 
and r]/X = 0.1 (solid curves). The HCM is a biaxial bianisotropic medium. 
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Figure 2: Real (left) and imaginary (right) parts of the HCM constitutive parameter £jjcm plotted 
against the eccentricity parameter p £ (0,1) for rj/X = (dashed curves), rj/\ = 0.05 (broken 
dashed curves) and rj/X = 0.1 (solid curves). The HCM is a biaxial bianisotropic medium. 
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Figure 3: Real (left) and imaginary (right) parts of the HCM constitutive parameters plotted 
against orientation angle 92 € (0,7r/2) for i]/\ = (dashed curves), r//A () = 0.05 (broken dashed 
curves) and ry/A = 0.1 (solid curves). The HCM is a biaxial bianisotropic medium. 
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Figure 4: Real (left) and imaginary (right) parts of the HCM constitutive parameter e x HCM plotted 
against relative correlation length k L 6 (0,0.2) for r]/L = (dashed curves), r//L = 0.5 (broken 
dashed curves) and rj/L = 0.95 (solid curves). The HCM is a biaxial bianisotropic medium. 
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Figure 5: Real (left) and imaginary (right) parts of the HCM constitutive parameters ^ X jf CM plotted 
against volume fraction f a G (0, 1) for 77/ A = (dashed curves), rj/X = 0.05 (broken dashed curves) 
and 7]/\ = 0.1 (solid curves). The HCM is a Faraday chiral medium. 
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Figure 6: Real (left) and imaginary (right) parts of the HCM constitutive parameter ^ CM plotted 
against the eccentricity parameter p £ (0,1) for rj/X = (dashed curves), rj/\ = 0.05 (broken 
dashed curves) and rj/\ = 0.1 (solid curves). The HCM is a Faraday chiral medium. 
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Figure 7: Real (left) and imaginary (right) parts of the HCM constitutive parameters €^f CM plotted 
against orientation angle ip € (0,7r/2) for r)/\ = (dashed curves), rj/\ = 0.05 (broken dashed 
curves) and r)/X = 0.1 (solid curves). The HCM is a Faraday chiral medium. 
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Figure 8: Real (left) and imaginary (right) parts of the HCM constitutive parameter ^ g HCM plotted 
against relative correlation length k L 6 (0,0.2) for r]/L = (dashed curves), r//L = 0.5 (broken 
dashed curves) and rj/L = 0.95 (solid curves). The HCM is a Faraday chiral medium. 
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